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Введение. В статье показано, что цифровой двойник (элек- 
тронный паспорт) станка с ЧПУ разрабатывается как ки- 
берфизическая система. 

Цель работы — создание нейросетевых моделей, опреде- 
ляющих функционирование станка с ЧПУ, его производи- 
тельность и динамическую устойчивость при резании. 
Материалы и методы. Рассматриваются вопросы созда- 
ния математических моделей процессов механической 
обработки с использованием системы сенсоров и про- 
мышленного интернета вещей. Оценены методы машин- 
ного обучения, подходящие для реализации названных 
задач. Предложена нейросетевая модель динамической 
устойчивости процесса резания, позволяющая оптимизи- 
ровать процесс механической обработки на этапе техноло- 
гической подготовки производства. На основе подходов 
нелинейной динамики реконструированы аттракторы ди- 
намической системы резания и определены их фракталь- 
ные размерности. Выбраны оптимальные характеристики 
оборудования по входным параметрам и отладке планиру- 
емого технологического процесса на основе цифровых 
ДВОЙНИКОВ. 

Результаты исследований. Использование методов ма- 
шинного обучения позволило создать и исследовать 
нейросетевые модели технологических систем обработки 
резанием и программное обеспечение для их реализации. 
Показана возможность применения деревьев решений для 
задачи диагностики и классификации неисправностей 
станков с ЧПУ. 

Обсуждение и заключения. В реальном производстве тех- 
нология цифровых двойников позволяет оптимизировать 
режимы обработки с учетом технического и динамическо- 
го состояния станков с ЧПУ. Это обеспечивает высоко- 
точную оценку производственных мощностей предприя- 
тия при составлении производственной программы. Кроме 
того, на основе интеллектуального анализа данных систе- 
мы распределенных сенсоров можно выявить неисправно- 
сти оборудования в режиме реального времени. 
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Введение. Предложенная в [1] единая базовая управляющая платформа служит основой при разработке 
нового поколения систем управления технологическим оборудованием, а также позволяет совершенствовать 
системы числового программного управления (ЧПУ) для цифровых производств. Такой платформой могут 
быть открытые системы ЧПУ станков с большими вычислительными ресурсами и высоким быстродействием 
обработки крупной базы данных (БД), встроенными модулями нейропроцессоров и модулями связи с промыш- 
ленным интернетом. В последнем случае речь идет о возможностях использования облачных технологий обра- 
ботки больших объемов данных на серверах предприятий (локальные сети) и провайдеров. Все это создаст ос- 
нову интеллектуального управления для широкой номенклатуры технологического оборудования с ЧПУ. 

С девяностых годов ХХ века [2] авторы представленного исследования разрабатывают программные 
продукты для интеллектуального управления технологическим оборудованием. В частности, были предложены 
критерии оценки динамической устойчивости процесса резания, основанные на методах нелинейной динамики 
и фрактального анализа сигналов виброакустической эмиссии (ВАЭ). Авторы создали единую платформу рас- 
ширения функциональных возможностей различного технологического оборудования с системами ЧПУ. Она 
реализуется путем встраивания в системы ЧПУ модулей высокопроизводительных вычислений и глубокого 
обучения искусственных нейронных сетей с использованием технологий СОПЛА пУ191а. В рамках данной рабо- 
ты изучается технологическое оборудование с ЧПУ, оснащенное датчиками и задействующее облачные техно- 
логии для сбора и обработки информации по разработанным методикам. Такие системы будут рассматриваться 
как киберфизические [3]. 

Предполагается, что разработанная авторами платформа станет основой цифровизации на всех уровнях 
предприятия. Для этого она должна не только анализировать данные оборудования, систем, устройств, но и 
использовать полученную таким образом информацию, чтобы сократить время вывода новых изделий на ры- 
нок, повысить гибкость производства, качество продукции и эффективность производственных процессов. 

Цифровой двойник (412 а1 11) — новое слово в моделировании оборудования, техпроцессов и плани- 
ровании производств [3]. Это совокупность математических моделей, достоверно описывающих процессы 
и взаимосвязи на отдельном объекте и в рамках целого производственного предприятия с использованием ана- 
лиза больших данных (61° даа) и машинного обучения. 

Лидер в области использования цифровых двойников — компания З1етепз [3]. Согласно ее определе- 
нию, цифровой двойник — это ансамбль математических моделей. Они характеризуют различные состояния 
оборудования, технологических и бизнес-процессов во времени, в соответствии с текущими производственны- 
ми условиями. Среди таких математических моделей особое место занимают нейросетевые схемы, т. е. 
нейросетевая модель процесса или изделия — это его цифровой двойник [2]. 

Основной технологической единицей цифрового производства является киберфизическая система 
(КФС) [1]. Ее высокие адаптивные и интеллектуальные возможности обусловлены следующими особенностя- 
ми: 

— ассоциативное восприятие информации, 
— постоянное обучение, 
— оценка текущего и прогнозирование будущего состояния. 

КФС способна анализировать многомерные данные, учитывая даже скрытые факторы реального произ- 
водства. Основываясь на этих данных, она может автономно решать задачи оптимизации и принимать правиль- 
ные решения. Поэтому КФС — главный элемент при создании цифрового двойника. Под цифровым двойником 
в данном случае подразумевается: 

— комплекс математических моделей, характеризующих различные состояния оборудования, технологических 
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— детализированные сборочные ЗО)-модели объектов, отражающие связи и взаимодействия между узлами. 

С этой точки зрения цифровой двойник можно рассматривать как цифровую идентичность КФС, элек- 
тронный паспорт, в котором фиксируются все данные об обрабатываемых материалах, выполненных техноло- 
гических операциях и испытаниях. 

В настоящее время цифровые двойники в основном создаются для коммерческого использования. Они 
успешно работают в нефтяной и газовой промышленности. При этом в литературе нет данных о двойниках 
оборудования механообрабатывающих производств на базе КФС. 

Цифровые двойники на базе КФС можно получить, используя: 

— традиционные аналитические подходы, основанные на математическом описании физических процессов; 
— современные статистические методы, в том числе машинное обучение. 

Материалы и методы. Методы машинного обучения, применяемые для построения статистических 
моделей, можно разбить на три группы: модели регрессионного анализа, модели классификации и модели вы- 
явления аномалий [4—6] (табл. 1). 

Таблица 1 
Основные методы машинного обучения 


Квантильная регрессия быстрого леса Самоорганизующиеся карты Кохонена 


Регрессия Пуассона Опе-у5-АП 


Выбор метода машинного обучения зависит от размера, качества и природы данных, а также от типа 





решаемых задач. Существующие методы требуют различных вычислительных мощностей и обладают различ- 
ной степенью точности. Как правило, их оценивают по возможности достичь точной аппроксимации данных и 
выявить границы в пространстве данных. Метод искусственных нейронных сетей (ИНС) является наиболее 
универсальным и точным, он позволяет работать с большим количеством данных и строить нелинейные зави- 
симости. Нейронные сети задействуют большое количество параметров настройки, что открывает возможности 
для создания высокоточных моделей процессов, работающих в режимах регрессионного анализа, классифика- 
ции и выявления аномалий. 

Основной метод анализа качества моделей — перекрестная проверка. Она позволяет оценить статисти- 
ческое качество исходных данных за счет построения и сравнения нескольких моделей, полученных на различ- 
ных обучающих и проверочных выборках. При построении моделей сложных объектов и систем требуется 
снижение размерности данных и устранение эффекта мультиколлинеарности переменных. Решение таких задач 
возможно за счет применения метода главных компонент, позволяющего представить многомерные данные в 
виде ограниченного числа составляющих (компонент). Такой обобщающий подход может применяться и для 
устранения переобучения моделей. Основные методы машинного обучения представлены в табл. 1. 

Для повышения качества модели используют бэггинг и бустинг. Суть данных алгоритмов заключается 
в построении не одной модели, а их комплекса, работающего над той же задачей. Результатом такой работы 
является некая интегральная оценка вероятности какого-либо события. Данная оценка может быть представле- 
на как синергетический эффект действия группы моделей, каждая из которых по отдельности работает неудо- 
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зующих различные комбинации методов машинного обучения и прошедших различные стадии верификации и 
улучшения. 

Итак, в общем случае цифровой двойник на базе КФС — это многофакторная модель оборудования [1], 
включающая ансамбль электронных (т. е. нейронносетевых) моделей. В данном случае определяющими будут 
следующие нейронные сети: динамической устойчивости процесса резания, сил резания и шероховатости обра- 
ботанной поверхности. Силы резания вызывают упругие отжатия в системе «инструмент — деталь», что обу- 
словливает погрешности формы и размеров. 

Для решения данной задачи был разработан комплекс статистических моделей, использующих методы 
машинного обучения. Полученные модели являются основой для цифрового двойника токарного станка с ЧПУ. 
Они позволяют решать задачи регрессионного анализа для прогнозирования динамики процесса резания при 
различных условиях обработки, шероховатости обработанной поверхности и сил резания. Кроме того, с помо- 
щью этих моделей решаются задачи классификации для оценки текущего состояния станка. Схема процесса 
создания цифрового двойника представлен на рис. 1. 


Создание цифрового двойника оборудования с использованием 
методов машинного обучения и облачных технологий 
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Рис. 1. Схема процесса создания цифрового двойника станка с ЧПУ 


Рассмотрим более подробно нейросетевую модель динамики процесса резания. Обучающая выборка 
была получена на основе телеметрических данных, собранных в процессе механической обработки с использо- 
ванием распределенной системы датчиков сигнала ВАЭ, динамометра и технологии промышленного интернета 
вещей (шиза ифегпей ог 1155, ПоТ). Количество датчиков, их тип и пространственная ориентация определя- 
лись в соответствии с компоновкой оборудования. Наличие в системе большого числа датчиков с различной 
пространственной ориентацией объясняется неоднородностью материалов и конструкций, особенностями рас- 
пространения сигналов, а также возможностью перестройки колебательной системы в процессе эксплуатации. 
Таким образом, использование гетерогенной системы сенсоров позволило получить наиболее полную динами- 
ческую картину процессов в п-мерном пространстве «состояние — время». В качестве протокола приема- 
передачи данных в рамках сети ПоТ был использован стандартный протокол ТСРЛР, а также текстовый формат 
обмена данными ]ЗОМ (ТауазсирЕ ОБ]ес{ МоаНоп). 

Для предварительной обработки сигналов применялись аналоговые и цифровые полосовые фильтры, а 
также вейвлет-фильтры. Использование вейвлет-фильтров позволило: 

— исключить влияние шумовой составляющей в сигналах ВАЭ; 
— выделить периодическую и хаотическую составляющие на основе энтропийных показателей. 

Предварительная обработка данных проводилась на «пограничном» устройстве (Едэе), которое вычис- 
ляет значения параметров сигнала ВАЭ и формирует пакеты данных для отправки в виртуальное облачное хра- 
нилище. Построение математических моделей на основе методов машинного обучения требует структурирова- 
ния и разметки данных. Однако учитывая особенности технологии работы ПоТ, построение реляционных БД не 
всегда возможно, поэтому был применен подход на основе технологий М№оЗОГ. Для этого хранение и обработка 
данных были реализованы на виртуальном сервере, имитирующем работу вычислительного кластера. На этом 


сервере была развернута специальная распределенная, масштабируемая файловая система. 
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В качестве основы использовались свободно распространяемые утилиты, библиотеки и фреймворки 
проекта Надоор (это одна из наиболее успешных и распространенных технологий больших данных). В частно- 
сти, программная оболочка Надоор Соттоп управляет распределенной файловой системой НПБЕЗ и базой дан- 
ных НВазе [4]. Для выполнения распределенных вычислений и обработки больших объемов данных использо- 
валась система МарКедисе, обеспечивающая автоматическое распараллеливание и распределение задач на кла- 
стере. Для ускорения обучения нейросетевых моделей были задействованы графические процессоры пУ1А1а 
СОВА. 

Результаты исследования. Для создания статистической модели динамического состояния процесса 
резания выбрана рекуррентная ИНС с сигмоидальной функцией активации нейронов. В качестве обучающей 
выборки использованы данные, полученные в процессе работы промышленного оборудования. При обучении 
использовался алгоритм обратного распространения ошибки. По окончании процесса обучения проводилась 
верификация полученной модели на основе значений средней квадратичной ошибки. На рис. 2, а видно моно- 
тонное убывание ошибки. При этом ее распределение является нормальным и сосредоточено вблизи нулевого 


значения, что свидетельствует о хорошем качестве полученной модели. 
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Рис. 2. Оценка качества полученной модели: диаграмма распределения ошибки обучения (а); зависимость 
среднеквадратичной ошибки от итерации обучения (6) 


Полученная ИНС (рис. 3) состоит из 17 входных нейронов, воспринимающих информацию об условиях 
обработки и предыдущем динамическом состоянии. 
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Рис. 3. Нейросетевая модель динамической устойчивости процесса резания 


Выходной слой нейронов отражает информацию о динамическом состоянии системы резания. Скрытый 


слой образован 38 нейронами. Данная ИНС способна моделировать динамику динамической системы резания 


для различных условий обработки в заданные моменты. 


Для оценки адекватности нейросетевой модели было проведено моделирование процесса механической 
обработки на токарном обрабатывающем центре с ЧПУ. В качестве обрабатываемой заготовки использовался 


вал из нержавеющей стали 12Х18НТОТ, закрепленный с двух сторон по центру, как показано на рис. 4, а. 
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Рис. 4. Результаты нейросетевого моделирования: схема обработки (а), моделируемый процесс резания (6), реальный 
процесс резания (в), процесс резания после оптимизации (г) 

Режимы определены в соответствии с производственными нормами для обработки деталей данного ти- 
па на токарном станке с ЧПУ. Симуляцию перемещения инструмента относительно заготовки обеспечила пе- 
ременная на входном слое ИНС, отвечающая за положение суппорта станка по оси Х. Шаг изменения данной 
переменной составлял 0,01 мм на расстоянии от 0 до 350 мм. Так было получено динамическое состояние си- 
стемы резания для 35 тыс. положений инструмента на оси 7. Зафиксированные динамические состояния отра- 
жены на соответствующей спектрограмме (рис. 4, 6): заметны две доминирующие частоты автоколебаний 
(3,5 кГци 6,5 кГц) и два фазовых перехода. 

Для подтверждения результатов моделирования проводился натурный эксперимент с аналогичными 
условиями обработки и регистрацией сигнала ВАЭ. По полученным экспериментальным данным была также 
построена спектрограмма процесса механической обработки (рис. 4, в). На ней тоже видны две доминирующие 
частоты в областях 3,5 кГц и 6,5 кГц и два фазовых перехода. В итоге погрешность моделирования на различ- 
ных участках обработки составила от 3 % до 7 %. Как видно из рис. 4, предложенные производственные нормы 
не позволяют получить изделие заданного качества. Это может быть причиной брака, а учитывая биения заго- 
товки, и поломки инструмента. 

Использование названных выше норм ведет к неадекватной оценке производственных мощностей 
предприятия и, соответственно, к некачественному формированию производственной программы. Устранение 
данного недостатка на этапе конструкторско-технологической подготовки чревато дополнительными времен- 
ными и финансовыми затратами. Поэтому выполняется корректировка режимов, заложенных в нормах произ- 
водства. Для этого рассчитываются оптимальные значения вектора технологических условий с использованием 
ИНС и метода градиентного спуска. Целевыми функциями оптимизации были значения амплитуды сигнала 
ВАЭ, фрактальной размерности и информационной энтропии, которые должны стремиться к минимуму. В ка- 
честве параметров оптимизации выбраны скорость обработки и значение фаски на передней поверхности. По 
завершении процесса оптимизации были получены новые режимы обработки: Г = 90 м/мин; 5 = 0,1 мм/об; 
у = 0,2 мм. 

Эти режимы оценивались также по результатам натурного эксперимента с регистрацией и обработкой 
сигнала ВАЭ. Как видно из рис. 4,г, оптимизация условий обработки позволила повысить динамическую 
устойчивость процесса резания за счет снижения амплитуды автоколебаний в 2 раза и фрактальной размерно- 
сти О{= 1,1. Таким образом, повышалось динамическое качество механической обработки без потери произво- 
дительности. 

Аналогичным путем были получены модели для прогнозирования значений шероховатости обработан- 
ной поверхности и сил при резании. Силы резания определяют точность обработки, вызывая упругое отжатие 
детали и инструмента, а, следовательно, и погрешность формы детали. Данные нейросетевые модели представ- 


ляют собой многослойные искусственные нейронные сети прямого распространения (рис. 5). 
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Рис. 5. Распределение значений амплитуды сигнала датчиков: А$_1 (а), Аз_2 (6), Аз 3 (в) 
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Нейросетевые модели в составе цифрового двойника применимы не только в реальных производствен- 
ных условиях, но и при выборе и обосновании закупки оборудования. На этом этапе возможен подбор опти- 
мальных характеристик оборудования и отладка планируемого технологического процесса на цифровых двой- 
никах. 

Для оценки текущего состояния оборудования была разработана модель классификации с использова- 
нием деревьев решений [5]. В качестве признаков, на основе которых принималось решение, были выбраны 
значения амплитуды сигнала ВАЭ в различные моменты времени с трех датчиков, расположенных на инстру- 
менте (Аз_1), у передней опоры шпинделя (Аз_З) и у заднего центра (Аз_2). В качестве классов использовалось 
три состояния: нормальное (№ шотта]), износ инструмента (Тоо] \меаг) и аномальное (АБпогта|). Под аномальным 
подразумевается любое состояние, при котором качество обработанной поверхности и работа оборудования 
будут неудовлетворительными. Критерием износа инструмента считалась фаска на задней поверхности 
й = 0,15 мм. После получения обучающей выборки выполнялся разведочный анализ данных. 

Из рис. 5 видно соответствие распределения значений амплитуды сигнала ВАЭ датчика Аз 1 трем 
группам состояний. Информация, полученная с этого датчика, хорошо дифференцирует нормальное состояние 
и износ инструмента. Однако значения, соответствующие аномальному состоянию, имеют пересечения с дан- 
ными из других групп, поэтому информация с датчика Аз$_1 не может быть использована с высокой точностью 
для классификации состояний. 

Рис. 5, б характеризует распределение значений амплитуды сигнала с датчика Аз_2. Полученные дан- 
ные имеют низкую информационную ценность из-за двух недостатков: 

— пересечения значений группы аномального состояния с другими классами; 
— слабая дифференциация нормального состояния и износа инструмента (об этом свидетельствует их близкое 
расположение). 

На рис. 5, в показано распределение значений амплитуды сигнала с датчика А$_3. Информация с этого 
датчика хорошо выделяет значения аномального состояния, при этом отсутствуют пресечения значений всех 
классов состояний. Таким образом, разведочный анализ данных позволяет сделать заключение о возможности 
получения статистической модели классификации состояния станка с ЧПУ. 

Для обучения модели был использован алгоритм САБКТ [3]. Он делит полученные данные на опреде- 
ленные подмножества на разных уровнях. В качестве критерия разделения используется целевая функция ми- 
нимизации информационной энтропии. Таким образом, каждый новый уровень имеет более определенное со- 
держание признаков. Совокупность всех полученных уровней представляется в виде иерархической древовид- 
ной структуры, состоящей из узлов и ветвей, при этом ветвление следует стандартным логическим правилам. 

Из всей выборки признаков выбирается такое значение условия деления для корневого узла, которое 
лучше минимизирует значения информационной энтропии. Затем от каждого полученного узла строятся ветви 
и происходит аналогичное деление на новые узлы. Подобное ветвление может продолжаться, пока не будут 
классифицированы все значения обучающей выборки. В общем случае такое дерево будет представлять слож- 
ную ветвистую структуру, идеально классифицирующую текущую выборку. Однако оно может показывать 
очень низкую точность при работе с новыми данными. Поэтому необходима оптимизация размера дерева. Она 
может проводиться с использованием перекрестной проверки, на основе которой строится несколько деревьев- 
претендентов, а в итоге выбирается дерево, показавшее лучший результат. Кроме того, выполняется отсечение 
ветвей, пока данная процедура не приведет к резкому возрастанию ошибки. 

На рис. 6 представлена полученная модель классификации состояния оборудования на основе показа- 


ний датчиков ВАЭ с использованием дерева принятия решений. 
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Рис. 6. Дерево решений для классификации состояния оборудования 


Первое выбранное условие — показания амплитуды сигнала ВАЭ с датчика Аз _3. Если амплитуда ВАЭ 
больше или равна 2,3, то состояние однозначно аномальное. Это суждение подтверждается разведочным анали- 
зом данных. В противном случае происходит следующее разбиение, и классификацию определяют данные с 
датчика Аз _1. Если амплитуда ВАЭ датчика Аз_1 меньше 2,1, состояние считается нормальным; если больше, 
то речь идет об износе инструмента. 

Как видно из рис. 6, принятие решений не учитывает показания датчика Аз_2. Низкая информатив- 
ность сигнала с этого датчика была выявлена и в ходе разведочного анализа. 

Высокая точность определения аномального режима работы на основе анализа сигнала датчика Аз$_3 
может быть связана с местом его установки (передняя опора шпинделя). Вероятнее всего, динамические про- 
цессы в шпиндельном узле обусловили это состояние, которое отражалось и на качестве механической обра- 
ботки. Учитывая нелинейные процессы рассеивания колебаний в упругой системе станка, данная информация 
доходила до других датчиков уже сильно искаженной. 

При этом датчик Аз _1, установленный на инструменте, позволяет выделять его износ с высокой точно- 
стью. Это можно объяснить доминирующей ролью резца в общей колебательной системе. 

Датчик Аз_2, установленный на заднем центре, с наименьшей эффективностью регистрирует информа- 
цию для классификации состояния. Это можно объяснить удаленностью Аз_2 от основных динамических про- 
цессов. 

Таким образом, описанный метод можно применить не только для задачи классификации состояния 
оборудования, но и для определения оптимального числа датчиков и мест их размещения. Однако стоит учиты- 
вать возможность перестройки динамических процессов в системе резания при изменении условий обработки и 
естественного износа деталей станка. 

Обсуждение и заключения. Использование цифровых двойников оборудования при планировании 
производства позволяет выявить узкие места в технологических операциях, повысить качество продукции, сни- 
зить риски поломки инструмента и нештатной работы оборудования. С помощью цифровых двойников можно 
оптимизировать режимы обработки с учетом технического и динамического состояния каждой производствен- 
ной единицы. Такой подход обеспечивает высокоточную оценку производственных мощностей предприятия 
при составлении производственной программы. Кроме того, на основе интеллектуального анализа данных си- 
стемы распределенных сенсоров в режиме реального времени выявляются неисправности оборудования. 
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